Inhomogeneous hard-sphere mixtures: Manifestations of structural crossover 
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We study various manifestations of structural crossover in the properties of a binary mixture 
of hard-spheres. For homogeneous mixtures that are sufficiently asymmetric, there is a crossover 
line in the phase diagram such that for thermodynamic states that lie on one side, the decay of 
the three bulk pair correlation functions is oscillatory with a common wavelength approximately 
equal to the diameter of the smaller spheres, and for states on the other side, the common wave- 
length is approximately the diameter of the bigger spheres. Using density functional theory we show 
that structural crossover manifests itself in the intermediate range decay of i) the one-body density 
profiles of a hard-sphere mixture adsorbed at planar walls, ii) the effective (depletion) potential 
between two big hard-spheres immersed in the same binary mixture and iii) the solvation force, 
or excess pressure, of the same mixture confined between two planar hard-walls. We determine 
exactly the structural crossover line for a one-dimensional binary mixture of hard-rods and present 
evidence, based on density functional theory calculations and Monte-Carlo simulations, for struc- 
tural crossover in homogeneous binary mixtures of hard-disks in two dimensions. By considering 
a multicomponent mixture of hard-spheres, with an appropriate bimodal distribution of diameters, 
we argue that structural crossover should still occur in the presence of polydispersity and that our 
results could be relevant to colloidal mixtures where correlation functions can be measured using 
real-space techniques. 

I. INTRODUCTION 

In 1969 Michael Fisher and Ben Widom pj published an elegant paper in which they conjectured that there should 
be a line in the phase diagram of a simple one-component fluid at which the character of the longest range decay of the 
total pair correlation function h(r) changes from monotonic to exponentially damped oscillatory. Their conjecture, 
which was based on exact results for certain one-dimensional models, pertains to fluids which exhibit liquid-gas 
phase separation. Although the Fisher-Widom (FW) line constitutes a sharp, well-defined crossover line, in say the 
temperature-density plane, it does not imply any accompanying non-analyticity in the free energy of the bulk fluid and 
therefore no phase transition. The significance of the FW line for the structure of wall-fluid and liquid-gas interfaces 
and for the asymptotic decay of the solvation force in a confined fluid was not appreciated until much later 0, . 
Indeed the first attempt at calculating the FW line using a realistic theory of liquids (for a square-well model) was 
not reported until 1993 2}. Subsequently there have been several investigations of FW crossover for various types 
of fluids using a variety of liquid state theories and the first computer simulation determination of the FW line (for 
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a truncated Lennard- Jones fluid) was reported in 2000 In the present paper we follow the spirit of FW in that 
we seek qualitative changes in the asymptotic decay of correlation functions corresponding to different regions of the 
phase diagram. Unlike FW, we focus on binary mixtures and investigate changes in the wavelength of the oscillations 
in various structural quantities; we are not concerned with crossover from monotonic to oscillatory decay rather with 
crossover from asymptotic oscillatory decay with one particular wavelength to decay with a different wavelength. 

In a recent paper denoted I, we showed that for sufficiently asymmetric binary mixtures of (additive) hard- 
spheres there is a sharp structural crossover line in the (rjb, r? s ) phase diagram at which the common wavelength of the 
longest range oscillations in the three bulk pair correlation functions hbb{r), h ss (r) and hb s (r) changes discontinuously 
from approximately Ub, the diameter of the bigger species, to <r s , the diameter of the smaller species, rjb and i] s denote 
the packing fractions of big and small spheres. We also showed, using density functional theory (DFT) and Monte 
Carlo simulation (MC), that structural crossover manifests itself at intermediate range, not just at longest range. 
It follows that crossover might be observed in real-space measurements of the pair correlation functions in colloidal 
mixtures. 

Here we investigate what repercussions this new crossover line might have for the structure of inhomogeneous 
hard-sphere mixtures. We also examine other properties of mixtures seeking further manifestations of structural 
crossover. 

Our paper is arranged as follows. In subsection III Al we summarise the main features of the general theory of 
the asymptotics of pair correlation functions in binary mixtures and outline the nature of structural crossover found 
in bulk fluids in paper I. Subsection III Bl presents DFT results for one-body density profiles, Pi(z), i — s,b, of the 
mixture at two types of planar wall. These demonstrate clearly that structural crossover is just as easy to observe at 
wall-fluid interfaces as in bulk. Indeed a single pole approximation for Pi(z) is just as effective as the corresponding 
approximation || for the pair correlation functions hij(r). In subsection IIII Al we consider the effective depletion 
W(r) potential between two big hard spheres immersed in the binary mixture. Using a particle insertion method, 
based on Widom's |6( potential distribution theorem, combined with DFT we show that the intermediate and long 
range behaviour of W(r) provides a clear signature of crossover, i.e. the wavelength of the oscillations in the depletion 
potential exhibits the same variation with the thermodynamic state point (rjb, Vs) &s is found in hij(r). An equivalent 
conclusion is reached in subsection llll Bl for the solvation force f s (L) which arises when the binary hard-sphere mixture 
is confined between two planar hard walls separated by a distance L. The oscillations in this thermodynamic quantity 
mimic closely those found in hij(r) and in the wall-fluid density profiles Pi(z). Throughout we focus on a mixture 
with size ratio q — u s /<Jb = 0.5. However, there is nothing special about that particular choice for q as we pointed 
out in paper I. 

In section iTVl we extend our analysis to the case of polydisperse hard-sphere mixtures, bearing in mind that poly- 
dispersity is always relevant in colloidal systems. We calculate, using DFT, the density profiles for a multicomponent 
mixture with a large, but finite, number of species adsorbed at a hard-wall. The size distributions are chosen to mimic 
a bimodal system with two maxima that are very weakly overlapping. In spite of the maxima being very broad we 
find a clear signature of crossover in the density profile of the species whose diameter is equal to a™ ax , the position 
of the maximum in the size distribution of big spheres. We return to homogeneous binary mixtures in section [V] 
and investigate structural crossover in low dimensional fluids. For one-dimensional hard-rod mixtures there are exact 
results for the pair direct correlation functions which allow us to determine exactly the relevant poles of the Fourier 
transforms of hij(r). We find crossover, demonstrating that this phenomenon arises in an exactly solvable model. 
The crossover line has a similar shape to that found for the hard-sphere mixture with the same size ratio. 

Two-dimensional hard-disk mixtures are of particular interest from an experimental point of view. Using video 
microscopy one can measure accurately correlation functions in colloidal suspensions that are confined by light fields 
to two dimensions. This technique might provide a means of testing the crossover scenario experimentally. We perform 
DFT calculations and Monte-Carlo simulations to verify that structural crossover also occurs for hard-disk mixtures 
in dimension d = 2. 

We conclude, in section l*VTl with a discussion of our results and their possible relevance for other mixtures and for 
experiment. 

II. ASYMPTOTIC DECAY OF CORRELATION FUNCTIONS 

In this section we recall some of the predictions of the general theory of the asymptotics of correlation functions 
laid out in Refs. 0, ; a summary of the main results is given in Ref. Q which adopts a density functional theory 
perspective from the outset. We consider a binary mixture of big (b) and small (s) particles with bulk (reservoir) 
densities pi and chemical potentials pi, i — s,b. Numerical results will be presented for hard-sphere mixtures. 
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A. Pair Correlations in Bulk Mixtures 



The nature of structural crossover of the total correlation functions hij(r) = gij(r) — 1, where gij{r) is the radial 
distribution function, was studied in detail for homogeneous binary hard-sphere mixtures in paper I. The analysis is 
based on calculating the leading-order poles of hij(k), the Fourier transform of hij(r). 

(2) 

In the bulk mixture hij(r) are related to the pair direct correlation functions c\j (r) via the Ornstein-Zernike (OZ) 
equations 



hi 



(r ia ) = i 2) (r 12 ) +Y,P»J d 3 r 3 c$(r 13 )h kj (r 32 ), (1) 



k—s,b 



where ry = [r^ — r 3 -|. 

If we suppose that the direct correlation functions are known the OZ equations can be solved easily in Fourier 
space, hij(k) can be written as 

hii (k) = (2) 
3K D(k) 1 ' 

where the denominator common to all three equations is given by 

D{k) = [1 - p s cf){k)][l - p b c$(k)} - PsPb c$(k) 2 . (3) 

c^\k) is the Fourier transform of c^'(r). 

From the formal solution one obtains an expression for the total correlation functions hij(r) in real space by 
performing the inverse Fourier transform via the residue theorem. We assume that the singularities of hij (k) are 
(simple) poles |2j. In order to determine the asymptotic behaviour we must consider the leading order pole (LOP) 
contribution to rhij(r) = (2tt)~ 1 R 1 ^ exp (ip n r), where p n are the aforementioned poles, given by the zeros of 

D(k). is the residue of khij(k) corresponding to the pole p n . Poles are either pure imaginary, p — iao, or occur as 
a conjugate complex pair p = ±ai + iao- The pole, or pair of poles, with the smallest imaginary part ao determines 
the slowest exponential decay and is termed the LOP. This pole determines the ultimate, r — > oo, decay of /iy(r). 
Note that in the case of complex poles all three total correlation functions show a common oscillatory asymptotic 
decay 

r hij(r) ~ Aij exp (— aor) cos (air — ©y), r — * oo. (4) 

The imaginary and real parts of the LOP determine the characteristic decay length Oq 1 and the wavelength of 
oscillations 2w/ai, respectively. These length scales are common for all pairs ij, whereas the phases 0,j and the 
amplitudes are dependent on the indices ij. As emphasised in I this observation is remarkable when one considers a 
binary system where the sizes of the two species are quite different. Clearly at high concentrations of the small particles 
one would expect the wavelength to be determined by the size of the small particles whereas at low concentrations 
one would expect that the wavelength of oscillations is set by the size of the big particles. Explicit calculations, based 
on DFT, show that the wavelength of the longest ranged oscillations 27r/ai changes discontinuously at some sharp 
crossover line in the phase diagram of binary hard-sphere mixtures [j| . In figure the crossover line is plotted for 
a hard-sphere mixture with size ratio q = a s /a b = 0.5. a s and a b denote the hard-sphere diameters and results are 
presented in terms of the packing fractions r\i = npiaf/6, i = s,b. 7ti denotes the pole with the longest wavelength 
27r/ai ~ (if,, whereas tt2 is that with the second longest wavelength In fa\ ~ <7{,/2. There are higher-order poles tt 3i 7T4, 
etc. 01 but these are not germane to figure The genesis of the crossover line is given in figure 2 of paper I where the 
trajectories of poles are plotted. Note that close to a crossover point the two poles 7Ti and ~ki have similar imaginary 
parts, ao and ao, and both contribute to the oscillatory decay of the total correlation functions at intermediate values 
of r, provided that the corresponding amplitudes, and Aij, are similar in size, i.e. near crossover we expect 

rhijir) ~ Aij exp (—aor) cos (air — Gy) + Ay exp (—aor) cos (air — 0y) (5) 

for large r. The first contribution corresponds to tt\ with a\ ~ 2w/ab and the second to 7T2 with cii ~ 27r/(<7<7&). 
Crossover occurs when ao = ah. 

In I we investigated the decay of hij (r) at intermediate range by performing DFT calculations of the density profiles 
of both species in the presence of a fixed (test) particle exerting an external potential on the other particles in the 
fluid. Using both the Rosenfeld Q and White Bear version [S| of fundamental measure DFT we found that the 
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FIG. 1: Crossover line for a binary hard-sphere mixture with size ratio q — 0.5. On each side of the line the asymptotic 
behaviour is dominated by different leading order poles (LOP) labelled m and 7T2. For small values of the packing fraction rj s , 
7Ti dominates and the wavelength of the oscillations in hij(r) is ~ cti,. On increasing r) s at fixed r/t the wavelength decreases 
gradually until crossover occurs. Then 7T2 dominates and the oscillations have a wavelength ~ 0.5<T(,. Generally at the crossover 
line the wavelength of the longest ranged oscillations changes discontinuously by a factor of about q. These results were obtained 
using the Percus-Yevick approximation for the direct correlation functions c'^'(fe) — see I. Note that the crossover line begins 
at the origin and we have truncated it at rji = 0.5. 



crossover behaviour predicted by the asymptotic (pole) analysis manifests itself at intermediate separations as well as 
at longest range, r — > oo. This is illustrated in figures 6 and 7 of I. Away from the crossover line the result from the 
LOP approximation provides an accurate fit to the DFT results except for small separations, r < 2ob, where many 
poles begin to contribute. For state points close to the crossover line the two-pole approximation J5| provides as good 
fit to the DFT results as is achieved by equation away from the line; the presence of two oscillatory wavelengths 
implied by equation JJjJ is manifest in the numerical DFT data. This observation will be important later when we 
consider manifestations of structural crossover for inhomogeneous fluids. 

We complete our summary of structural crossover in homogeneous binary hard-sphere mixtures by presenting, in 
figureEl results for the real part, ai, of the leading order poles calculated for a full range of size ratios q. As in I, the 
poles were determined by calculating the zeros of D(k) in equation @ using the Percus-Yevick results for the pair 
direct correlation functions. For the more symmetric mixtures, q 0.7, the LOP is 7Ti for all packing fractions so 
there is a continuous variation of a\ in this regime; there is no crossover. For more asymmetric mixtures there are 
discontinuities in a\ and separate branches appear in figure [21 For example, for q = 0.5 there is a gap between the 
branches corresponding to 7Ti and 7r2, and crossover corresponds to jumping between these branches as the packing 
fractions are varied. As q is reduced further additional branches appear and for q = 0.35 crossover can occur between 
7Ti and 7T2 and between 7r 2 and 1T3 or directly between tti and 773. Similarly for q = 0.2 crossover can occur between 
7Tx and 7T4, 7T4 and 7r5 or directly between 711 and 7T5- 

The origin of the crossover lines shown in figure 5 of I can be understood from this 'master' plot of the poles. Note 
that for small size ratios some of the packing fractions considered in the plot correspond to state points that lie within 
the fluid-solid coexistence region — see I. 



B. Density Profiles at a Planar Wall 

Earlier treatments 0, 0, °f the asymptotics of correlation functions emphasised that, for fluids with short- 
ranged interatomic potentials, the pole structure given by the zeros of D{k) in equation J3J determines not only the 
decay of hij (r) in a bulk mixture but also the decay of the one body density profiles at wall-fluid and liquid- vapour 
interfaces. Given this observation we should expect to find manifestations of structural crossover in the properties of 
inhomogeneous mixtures and we seek these in the present subsection. 

Perhaps the most direct way of analysing the asymptotic decay, z — > 00, of the density profiles, p W i(z), of a binary 
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FIG. 2: Real part a± of the leading order poles calculated for fixed size ratios q and various combinations of packing fractions: 
generally r\ — n 3 + rjb increases from left to right up to a maximum of 0.5. For the less symmetric mixtures with q < 0.65 
separate branches tti, tt2, etc. arise and there are no leading order poles with values of ai lying in the gaps. Structural crossover 
corresponds to a 'jump' from a value of ai associated with one branch of poles to a value associated with another — see text. 



mixture adsorbed at a planar wall is to start from the wall-particle OZ equations, i.e. we consider a ternary mixture 
with species s,b plus a third species whose density ps — > 0. The third component, whose diameter is 03, can be 
considered as a wall that exerts an external potential on the mixture of s and b particles. The resulting OZ equations 
can be expressed as 

h wi {k) = ^777- , (6) 
D(k) 

i = s,b, where h W i(k) denotes the Fourier transform of h W i(r), the wall-particle total correlation function and the 
third species is now labelled w (wall). Note that the form of © is valid for any value of 03 and the repercussions for 
the density profiles hold for a spherical wall, 03 < 00, for which the wall-particle total correlation function h W i{r) has 
spherical symmetry, as well as for a planar one, 173 — > 00, for which we have planar symmetry for the total correlation 
functions, h W i(z). z is then the distance from the planar wall and h W i(k) is replaced by h W i{k) 1 the one- dimensional 
Fourier transform. The denominator D{k) is common to both species i = s and b particles and is identical to that 
appearing in equation it depends upon the pair direct correlation functions of the bulk mixture of s and b whose 
densities p s and pt, are fixed by the reservoir far from the wall. In the planar case the numerator N wi (k) can be 
written in terms of the one-dimensional Fourier transforms of the wall-particle direct correlation functions c W i(z) 
which depend, of course, upon the wall fluid potentials V^(z). Provided the latter decay faster than correlation 
functions in the bulk fluid the asymptotic decay of h W i(z) will be determined be the LOPs of D(k). Since bulk 
correlations decay as exp (— a r), where a is the imaginary part of the LOP, this requires V^(z) to decay faster than 
exp(— arj-z). In the calculations to be described below we shall consider only hard walls or walls whose attractive 
potential is of finite support. 

The planar wall-particle density profile is given by p W i{z) = Pi[h W i(z) + 1], i = s, b, and from equation JSJ) it follows 
that the asymptotic decay is described by 

pwi{z) — Pi ~ Bf exp (— aoz) cos (a±z — Qf), z — * 00 (7) 

where pi denotes the bulk density of species i. Although the decay length a^ 1 and wavelength 2it ja\ are the same for 
both species, the amplitude Bf and phase Qf depend on the species label and on the choice of wall-fluid potentials. 
Clearly equation (J7J is the analogue of equation Note that the additional factor 1/r in the latter arises from 
taking a three-dimensional rather than a one-dimensional Fourier transform. 

Equation Q implies that if structural crossover occurs in the correlation functions of a bulk mixture one should 
observe equivalent behaviour in the one-body density profiles of the same mixture, near the same crossover state 
point, adsorbed at a planar wall. We examine this prediction below. 
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FIG. 3: Logarithm of density profiles p w b(z) of the big species b of a binary hard-sphere mixture, with size ratio q — 0.5, 
adsorbed at a planar hard wall. The bulk packing fraction r/b is fixed at 0.1. The solid lines refer to DFT calculations and 
the dashed lines result from leading order asymptotics — see text. For a low packing fraction of the small species n s = 0.1 
[curve (a)] the oscillations have a wavelength of about ab, the diameter of the big spheres whereas for a high packing fraction 
rjs = 0.15 [curve (c)] the wavelength is about a s , the diameter of the small spheres. In the vicinity of the crossover, curve (b) 
with rjs = 0.125, interference effects are clearly evident. In (a) and (c) a contribution from a single (dominant) pole, equation 
J7j, is used in the asymptotic expression whereas in curve (b) the contributions of the two poles are included, equation 
The results in (b) and (c) are shifted vertically for clarity of display. 



There are alternative ways of deriving equation based on i) a DFT approach [7j, or ii) analysis of the exact 
integral equations for the density profiles in planar geometry |Toj . 

In our numerical calculations we consider the binary hard-sphere mixture for which we reported the crossover line 
in figure^] namely with a size ratio q — 0.5. We determined the density profiles p W i(z) using Rosenfeld's fundamental 
measures DFT |8(, treating the planar wall as an external potential for the fluid mixture. Two types of wall- fluid 
potential were investigated. The first corresponds to a planar hard-wall: 

V ^ = { 0,' otherwise, ^ 

for i = s, b. Here z denotes the centre of a hard sphere, Rb = erf,/2 and R s = a s /2. In figure|21we show density profiles 
for the big spheres, p w b{z), for a fixed bulk (reservoir) packing fraction rjb — 0.1 and three values of n s , the bulk 
packing fraction of the smaller spheres. For n s = 0.1, curve (a), the oscillations have a wavelength of approximately 
ab whereas for rj s — 0.15, curve (c), the wavelength is approximately a s — qob- In curve (b), r\ a = 0.125, the density 
profile exhibits interference effects and two wavelengths are clearly evident. This value of n s lies close to the crossover 
value r\* = 0.126 for this choice of r\b — see figure^ Structural crossover manifests itself in a very direct fashion. The 
other striking feature of these results is that the oscillatory behaviour predicted by leading order asymptotics sets in 
at intermediate distances, i.e. for z > 2.5(7;,. This is illustrated by plotting alongside the DFT results the prediction 
from equation Q, with ao and a\ calculated from the zeros of D(k) in equation @ corresponding to the appropriate 
bulk state point (rjb,r] s ). As in subsection III Al the input to D(k) are the pair direct correlation functions cffi(k) 
obtained from Percus-Yevick theory. This is completely consistent with employing the Rosenfeld DFT since in bulk 
the latter generates, via functional differentiation, the Percus-Yevick c[^\r) 0,0- 

The dashed lines in curves (a) and (c) of figure [3] are results based on the single pole approximation, equation (7J, 
with the amplitude and phase fitted to the DFT results at intermediate z. The agreement between the two 
sets of results for p w b{z) is excellent, apart from the region very close to the wall where many poles contribute. For 
curve (b), the state point close to crossover, we utilise a two pole approximation, equivalent to equation (JSJ): 

Pwi(z) - Pi ~ Bf exp(— aoz)cos(a 1 z - Qf) 

+Bf exp (— a z) cos {a\z — <9™), z — > oo (9) 
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i = s,b. Once again ao, ao, a\ and Si are obtained from the zeros of D(k) at the given state point and now two 
amplitudes B™ and B™ and two phases 6™ and 6™ are fitted to the DFT data at intermediate z. The two-pole fit 
clearly provides an excellent description in the crossover regime. Overall the quality of the fits achieved by leading 
order asymptotics is no worse than is achieved by the corresponding fits to the test particle DFT results for hij (r) 
obtained for size ratio q = 0.3 — see figure 6 of paper I. The implication of this observation is that the crossover 
features found in the DFT results should be clearly visible in computer simulations of hard wall-fluid density profiles 
which, owing to statistical considerations, are limited to small and intermediate z. 

The second type of planar wall that we consider is described by the wall-fluid potentials 

!oo, z < Ri 

-e+^(z-R i )-^(z-R i ) 2 , RiKzKX + Ri (10) 
0, otherwise, 

where (i = (keT)- 1 and e > is a dimensionless constant describing the depth of the attractive well, i.e. PVp AW (Ri) = 
—£. The range of the potential is given by the length scale A and the form is chosen so that the potential and its first 
derivative are zero at z = A + Ri. In our calculations we set A = 1.25(76 and e = 1. 

Including such an attractive component to the wall-fluid potentials will modify the density profiles in the vicinity 
of the wall but should not alter the fundamental character of the long ranged decay; leading order asymptotics should 
be just as accurate as in the case of the purely hard wall. Of course, the amplitude and phase in equation J7J will 
be altered. In figure 01 we show the density profiles p W i(z), i = s,b, of the binary hard-sphere mixture with q = 0.5 
calculated for the bulk state point r/b = 0.1, rj s = 0.1. Results for the hard- wall, described by equation JHJ, are 
compared with those for the attractive wall, equation 1)10(1. It is energetically favourable to have a large number of 
particles in the attractive well. Since small spheres can pack better close to the wall than the larger ones it follows 
that the contact value, p s (a s /2), is considerably larger for small particles at the attractive wall than is the case for 
the same species at the hard- wall. Figure 0] demonstrates that adding an attractive piece to the hard- wall potential 
has a pronounced effect on the form of the density profiles of both species. From figure 0] one cannot observe any 
commonality of structure between the profiles of big and small species or between profiles corresponding to the two 
choices of wall-fluid potentials. Rather it is necessary to study the density profiles at intermediate and long range in 
order to observe such common features. In figure we plot the logarithm of the profile of the big spheres for z up 
to 10(7;,, the limit of our numerical accuracy. As was the case for the hard-wall, the asymptotic behaviour appears 
to set in at intermediate distances, z > 3crf,. Away from crossover, e.g. in curve (a) with rj s = 0.075 and curve (d) 
with T] s = 0.15 the results for the two types of wall-fluid potential are very close. For both types the wavelength of 
oscillations is close to in (a) and close to a s in (d). Clearly the asymptotic decay remains dictated by a single 
pole contribution when the wall exerts some finite ranged attraction on the fluid. For values of rj s lying closer to 
crossover the situation is somewhat more complex. Curve (b) in figure |5j for rj s =0.1, shows that the attractive wall 
potential influences strongly the form of the density profiles out to distances z ~ 7<Jb and the asymptotic behaviour, 
characterised by the wavelength cr&, does not set in until larger z. For rj s = 0.125, curve (c), there are clear differences 
between the results for the two different wall potentials extending out to z ~ 10<r&. It would appear that for the 
attractive wall the intermediate range behaviour of the profile is already dominated by the length scale a s , although 
crossover does not occur until the slightly larger value rj* — 0.126. Clearly the amplitudes entering the two-pole 
approximation |J5J depend sensitively on the choice of wall-fluid potential. Presumably for the attractive wall the 
amplitude B™ corresponding to the pole with a% — 2tt/<t s is much larger than B™, that for the pole with a\ = 27r/<7b 
whereas the imaginary part of the pole ao is slightly less than Oo- 

We conclude that crossover in the wavelength of oscillations for density profiles is not restricted to purely repulsive 
walls but can also be observed when the wall exerts attractive forces on the fluid. 

Note that although we have focused on the profile of the big species, our results for the small species show equivalent 
features — as is implied by equations (7J and JjJJ. We have also investigated density profiles for other size ratios q 
and we find similar manifestations of crossover near the crossover point appropriate to the given value of q. Overall 
the behaviour is equivalent to that found for hbb(f) in paper I. 

III. DEPLETION POTENTIALS AND SOLVATION FORCES 

In this section we turn attention to the asymptotic decay of two, closely related, properties of liquids. These are 
i) the depletion or solvent mediated potential, W(r), which arises between two big particles immersed in a solvent 
and ii) the solvation force, f s (L), which is the excess pressure arising from confining the solvent between two planar 
walls. We present evidence for structural crossover in W(r), at large particle separations r, and in f s (L), at large wall 
separations L, when the solvent is the binary hard-sphere mixture described in earlier sections. 
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FIG. 4: DFT results for the density profiles p w i(z) of the components i = b, s of a hard-sphere mixture adsorbed at the planar 
attractive wall, described by equation HUH , compared with the profiles at the planar hard wall, equation (HJ, for the case 
rit — 0.1, rj s = 0.1 and q = 0.5. The additional short-ranged attraction of the wall ensures that the small particles (s) exhibit 
a higher density at contact, p ws (a s /2), and for z close to the wall than they do at the hard wall. For the short distance scale 
shown here the density profiles clearly reflect the influence of the attraction. 
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FIG. 5: As in figure^Kg = 0.5, 776 = 0.1) but now we plot the logarithm of the density profile p w b(z) of the big species b over 
a much larger range of z. (a) rj s = 0.075, (b) rj s = 0.1, (c) rj a = 0.125 (close to the crossover value) and (d) rj s = 0.15. The full 
lines refer to results for the attractive wall, equation l|10|l . and the dashed lines to those for the hard- wall, equation |jSJ. In (a) 
and (d) the pair of profiles are very similar for z > 3<7(,, whereas in (b) the asymptotic behaviour for the attractive wall does 
not appear to set in until z > 7cr{,. In (c) it is difficult to resolve two distinct wavelengths for the attractive wall. The curves 
are shifted vertically for clarity of display. 



A. Depletion or Solvent Mediated Potential 



Consider two big particles of species 3 immersed in a solvent. Upon integrating out the degrees of freedom of the 
solvent atoms or molecules one obtains an effective two-body potential 

4>g (r) =0 33 (r)+W(r), (11) 
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where ^33 (r) is the bare interaction between the solute particles and W(r) is the solvent mediated potential arising 
from solute-solvent and solvent-solvent interactions. W(r) depends on the thermodynamic state of the solvent. A 
well-known example is the case of two hard-spheres, modelling colloids, immersed in a sea of ideal, non-interacting 
particles, modelling 'ideal polymer'. The exclusion of the 'polymer' from the hard spheres gives rise to a depletion 
zone between the spheres and the resulting W(r) is the celebrated depletion potential introduced by Asakura and 
Oosawa ^10] and by Vrij 0|. For the ideal solvent W(r) is attractive and monotonically increasing with r up to 
a certain finite separation beyond which it vanishes. When the solvent is a hard-sphere fluid and the solute particles 
are again (big) hard-spheres, W(r) exhibits depletion driven attraction at small separations and oscillatory decay at 
intermediate and long range. The oscillations arise from the packing of the small spheres that constitute the solvent 
and their wavelength is determined by the diameter of the small species, e.g. Ref. [l4j . 

Here we investigate the depletion potential when the solvent is the binary mixture of s and b hard-spheres. On 
intuitive grounds one might expect the oscillatory structure in W(r) to reflect the form of correlations in the bulk 
binary mixture and, therefore, to show crossover behaviour. In order to make this explicit we generalise the argument 
given in Ref. for a one-component solvent to a two-component solvent. For a ternary mixture of species s,b and 
3 the Fourier transform of the total correlation function h^(k) is given by an expression equivalent to equation JSJ. 
In the limit where the density of species 3 vanishes, p% — ► 0, the denominator should reduce to that in equation 
as was the case for wall- particle correlations in equation © , so the asymptotic decay of /133 (r) will be determined by 
the zeros of D(k) at the appropriate state point. In the same limit we have 

M»0 + 1 = 933(0 = exp [-/3<f$ (r)]. (12) 
It follows that for large r, where we can linearise the exponential, 

A33 

- f3W(r) ~ h 33 (r) ~ exp (-a r) cos (air - 633), r -> 00. (13) 

r 

Note that we assume the bare potential ^33 (r) is of shorter range than W(r). The same LOP that determines the 
decay of hij(r), i = s,b, in the bulk mixture must determine the decay of the depletion potential. Of course, the 
amplitude A33 and phase ©33 depend upon the solute-solvent interactions. 

In order to determine whether structural crossover can be observed in the depletion potential we performed calcula- 
tions of W(r) using the particle insertion method employed previously in DFT studies for one-component hard-sphere 
solvents |l4j . We first fix a single sphere of diameter 03 at the origin and make this an external potential for the fluid 
mixture of s and b hard-spheres. As discussed earlier, the density profiles of both s and b particles, obtained using 
the Roscnfcld DFT, display structural crossover as r — > 00. In a second step, a second particle of species 3 is inserted 
into the inhomogeneous mixture of s and b and input from the three component DFT is used to calculate W(r) |14|. 
The results shown in figure |H1 refer to the hard-sphere mixture with q = 0.5 and r/t, = 0.1 considered earlier. We set 
03 = 5er s . At short distances from contact (h = in figure El where h = r — 173 is the separation between the surfaces 
of the two big spheres of species 3) the depletion potential is attractive whereas oscillatory structure develops for 
larger distances. Results are presented for the same three choices, rj s =0.1, r/ s — 0.125 and r) s = 0.15, employed in 
figure^ From figure[|Ja) we note that the width of the attractive depletion well seems to be determined by the size of 
the majority species in the mixture. Thus, for rj s = 0.1 the big spheres dominate and the width of the attractive well 
is about 0.4<7b. On the other hand for r] s = 0.15, where the smaller species dominates, the attractive well is deeper 
and narrower and the depletion potential exhibits a pronounced potential maximum near h = OAab- 

In figure EJb) we show the same data plotted as In |/3VF(ft)| in order to expose the nature of the oscillatory decay. 
For rj s = 0.1 the wavelength is about Cb whereas for r/ s — 0.15 this is about a s — crf,/2. Interference is clearly visible 
for r) s — 0.125, very close to the crossover value. The sequence of curves is very similar to that displayed in figure |3| 
for the density profiles at a hard planar wall. We confirmed that the crossover mechanism is identical by plotting (as 
dashed lines) alongside the DFT results those obtained from the single pole approximation (|13fl in which ao and a\ 
refer to the LOP and are equal to the values employed in figure|3] For rj s — 0.1 and 0.15, i.e. curves (i) and (iii), the 
single pole approximation gives an excellent account of the DFT results for h > 3.5ab. Close to the crossover point, 
curve (ii) with r] s = 0.125, two LOPs are required, i.e. we employ the analogue of equation JSJ). As in figure|2| we fit 
amplitudes and phases at intermediate separations. By comparing figure E^b) with figure [21 we see that crossover is 
just as evident in the depletion potential as in the density profiles at a planar wall. Note that numerical noise sets 
into our results for h > 8<7b. 
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FIG. 6: Depletion potential W(h) between two hard-spheres of diameter 03 = 5a 3 immersed in a binary mixture of hard spheres 
with diameters a s and a^, = 2<r s , i.e. q — 0.5. h = r — oz is the separation between the surfaces of the two big spheres of species 
3 at distance r. The packing fraction of species b is rjb — 0.1. DFT results are shown in (a) for three different values of r) s . 
Increasing rj a reduces the width of the attractive well and increases its depth at contact h = 0. In (b) we plot the logarithm 
of over a larger range of separations. For r) s = 0.1 [curve (i)] the oscillations have a wavelength of about (Jb whereas 

for 778 = 0.15 [curve (iii)] the wavelength is about o 3 . In curve (ii) r\ a = 0.125, very close to crossover, interference effects are 
evident. The solid lines in (b) are the DFT results and the dashed lines result from leading order asymptotics — see text. The 
results (ii) and (iii) are shifted vertically for clarity of display. These results demonstrate the same crossover behaviour as the 
density profiles in figure [3] 



B. Solvation Force 

The solvation force is defined 0, 0] for a planar slit of width L as 

where f2 is the grand potential of the fluid mixture and p is the pressure of a homogeneous mixture with the (reservoir) 
chemical potentials {/ij}, i — s,b, and temperature T. We consider the limit where the area of the walls A — > 00. 
f s (L) is then the excess pressure resulting from confining the mixture in the slit; in the limit L — > 00, f s (L) — » 0. 
Here we consider only fluid-fluid and wall-fluid contributions to the solvation force, not the direct interaction between 
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the the walls ^(|. In the asymptotic regime, L — > oo, one might expect the decay of f s (L) to be determined by 
the properties of the bulk mixture described by the given {pi}, T. However, it is not immediately obvious that the 
thermodynamic quantity f s (L) should decay in the same fashion as the bulk pair correlation functions hij(r), i.e. with 
the same exponential decay length a^" 1 and wavelength 1%ja\. 

There are several ways of establishing this result. The most intuitive argument is to consider the ternary mixture 
described in the previous subsection, with p-$ — > and, in addition, to allow the diameter of species 3 to become 
macroscopically large. The effective potential between the macroscopic particles will still decay as in equation l|13fl . 
Attard and co-workers 17] placed this argument in the framework of the wall-particle OZ equations (see subsection 
III Bp and showed that the decay of the excess interaction free energy per unit area should be governed by the zeros 
of D(k), provided that the interatomic forces are sufficiently short-ranged. Such a result is also implicit in earlier 
integral equation studies by D. Henderson and co-workers, e.g. Ref. [Tsl ]. 

In the particular case of perfectly hard-walls, for which equation (|14fl reduces to the exact expression 

pf„{L) = iP*A0 + ) - Pi,oc(0 + )], (15) 

i— 6,s 

where Pi,L(0 + ) denotes the one-body density of species i at contact with one wall of the slit of width L and Pi tOO (0 + ) 
is the same quantity for infinite wall separations, there are additional arguments |2| to support the contention that 

f s (L) ~ Fexp (— aoL) cos {a\L — ©f)> L — > oo (16) 

where F and <3p are the amplitude and phase, respectively. Both quantities depend on the wall-fluid potentials. 

In the following we study the binary-hard sphere mixture confined by two planar hard walls described by the 
potential 

rri , s _ / oo, z < Ri and L - R t < z , . 

Valid*)- j o, R l <z<L-R l , [U > 

which is the generalisation of equation (jHJ. In our calculations we use the external potential equation (|17|) as input 
into the Rosenfeld fundamental measures DFT. We obtain the contact values Pll(Q + ) from the density profiles Pi,h{z) 
resulting from minimising the grand potential functional f2[{pi(z)}]. Note that the profiles obtained from the DFT 
satisfy the contact value theorem, i.e. the solvation force obtained from equation (|15|l is consistent with that obtained 
from differentiating w.r.t. L the equilibrium grand potential f2 resulting from minimising the functional. Considerable 
computational effort is required to produce sufficient values of the solvation force that we can have smooth curves as a 
function of the separation L. In order to investigate crossover behaviour we focused on the same mixture with q = 0.5 
and reservoir packing fraction = 0.1 and varied rj s over a similar range to that considered in earlier sections. 

In figure [7| we plot the logarithm of \(3aff s (L)\ for three values of rj s . For a small packing fraction, curve (a) 
T] s = 0.07, the oscillations in f s (L) have a wavelength of about df, for L > 3.5<7b. At a high packing fraction, curve (c) 
7] s = 0.15, the wavelength of oscillations is given by the diameter of the small spheres a s for L > 2ov At intermediate 
values of rj s we observe interference effects in the oscillatory decay. In curve (b), where rj s — 0.11, the oscillations 
at small L have a wavelength of about a s whereas at intermediate L interference effects resulting from two distinct 
contributions are apparent. Finally for large L, the oscillations have a wavelength of about o\>. Clearly confinement 
results in the small spheres imposing their length scale on the oscillations in f s (L) at small L and restricts interference 
effects to larger L. Recall that the crossover value for this choice of r)b is rj* = 0.126. Comparing the results in figure 
0with the density profiles p w b(z) for a single planar hard-wall in figure [3] we find close similarities. In particular, 
curves (b) and (c) in figure [7| provide compelling evidence that crossover has occurred in the form of oscillatory decay 
of f s (L) for packing fractions between rj s = 0.11 and t] s = 0.15. Note that there is some noise in the latter results for 
L > 9cf,, the limit of our numerical accuracy. 

It is important to recognise that there is no jump in the solvation force, or its derivatives at the crossover line. Near 
crossover two LOPs contribute to f s (L), in the manner of equations JSJ and JJjJ, and the real and imaginary parts of 
both poles vary smoothly with state point. 



IV. POLYDISPERSE HARD-SPHERE MIXTURES 



From an experimental point of view it is important to keep in mind that systems of colloidal particles with identical 
size (or shape) cannot be produced. It is not possible to produce a perfect binary mixture; there will always be 
a distribution of size ratios in even the best bimodal situation. Here we examine the influence of polydispersity 
on the structural crossover which was clearly manifest for a mixture of hard-spheres with a unique size ratio. The 
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FIG. 7: Logarithm of the solvation force f s (L) calculated from DFT for a binary hard-sphere mixture, with size ratio q = 0.5 
and 77b = 0.1, confined between two planar hard walls. Upon increasing the packing fraction of the small particles from rj a = 0.07 
[curve (a)] to r/ s = 0.15 [curve (c)] we find a crossover in the wavelength of oscillations. Interference effects can be observed at 
intermediate separations L for rj s = 0.11 [curve (b)]. Curves (b) and (c) have been shifted vertically for clarity of display. 



Rosenfcld DFT can treat mixtures with an arbitrary number of different components (species). However, we restrict 
consideration to a finite number of different particle sizes and replace each of the particle types big, 6, and small, s, 
by a 21-component hard-sphere mixture. 

The size distributions are chosen to be binomial and are shown in figure |H] In the case we consider the ratio of 
diameters ( j™ al /(7™ al i corresponding to the positions of the two maxima, is 0.5. The sum of the packing fractions, 
J2i Vb,i of the big particles is fixed at 0.1 while that of the smaller particles takes values J2i Vs,i — 0.1, 0.15, and 0.2 
in order to compare with the results from earlier sections. We treat rather broad distributions in order to capture 
effects which might arise in experimental situations. The width of both distributions is 0.6 and there is a very small 
overlap — see figure |S1 Figure shows the density profiles p w b(z) of the species with diameter a™ ax for polydisperse 
hard-sphere mixtures adsorbed at a planar hard wall. The Rosenfeld grand-potential functional is minimised and the 
resulting Euler-Lagrange equations for the profiles of the multi-component mixture are solved imposing the boundary 
conditions that the reservoir packing fractions are those given in figure [H] The results for p w b(z) are very similar to 
those in figure|3|for the binary hard-sphere mixture with size ratio q — 0.5. In curve (a), 7} S) j = 0.1, the wavelength 
of the oscillatory decay is characterised by the big spheres whose reservoir packing fraction rjb.i is a maximum. The 
small particles determine the wavelength of the oscillations in curve (c), ^iVs.i — 0.2. Curve (b), J^iVs.i — 0.15, 
clearly displays interference effects. If we take the presence of the latter to signal crossover then this is occurring 
for a larger (total) small sphere packing fraction than is the case for the binary mixture, where crossover occurs at 
rf s = 0.126. This observation is consistent with results from paper I where we found that for rjb = 0.1 the crossover 
value of t7s increases as q increases or decreases (see figure 5 of I), i.e. the minimum value of the crossover packing 
fraction occurs near q = 0.5. In the polydisperse mixture one is sampling component pairs with size ratios different 
from q = 0.5. 

It is striking that for such broad distributions there is still a clear signature of crossover in the density profile and 
we shall return to this point in the discussion. 



V. STRUCTURAL CROSSOVER IN ONE AND TWO-DIMENSIONAL HOMOGENEOUS FLUIDS 

All the results that we have presented so far refer to three-dimensional hard-sphere mixtures. Moreover these have 
been obtained by approximate methods, i.e. via DFT or integral-equation approaches. In this section we enquire 
whether structural crossover manifests itself in lower dimensional fluid mixtures. There are two distinct reasons for 

(21 

considering such systems. In one-dimension the pair direct correlation functions c i} ■ (r) can be calculated exactly for a 

homogeneous mixture of hard-rods of different lengths which implies that the poles of hij (k) can be determined exactly. 
The reason why two-dimensional systems are of interest is that real-space experimental techniques for colloidal fluids 
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FIG. 8: Distributions of the components in a bulk polydisperse hard-sphere mixture plotted as the ratio of the hard-sphere 
diameter a to a™ ax , the diameter of the component with maximum packing fraction in the distribution of the bigger spheres. 
Each binomial distribution consists of 21 components and the sum of packing fractions of the bigger components V_\ T]b,i = 0.1. 
The sum of the smaller components J^rfoi takes values 0.1, 0.15, or 0.2. Note that the ratio a™ ax /a™ ax = 0.5; 
is the diameter of the component with maximum packing fraction in the distribution of the smaller spheres. 



often restrict the determination of the pair correlation functions to two dimensions [2d I21L |2 



A. One-Dimensional Mixtures of Hard-Rods 

One-dimensional hard-rod fluids occupy a special position in the theory of liquids since many of their equilibrium 
properties can be calculated exactly for both homogeneous and inhomogeneous cases. In particular the exact intrinsic 
free energy functional is known for both the one component fluid and mixtures of hard-rods |24j . For our present 
purposes it is sufficient to note that Rosenfeld's fundamental measure DFT reproduces the exact functional when 
applied to the one-dimensional inhomogeneous hard-rod system [8|, l25j . The pair direct correlation functions are 
given by an expression equivalent to equation (11) in I but with different weight functions and a different excess free 
energy density /3<E>. Their (one dimensional) Fourier transforms c^\k) are given in Ref. [25| . 

We used these results to obtain the poles of hij (k) . These are given by the zeros of an expression equivalent to that 
for D(k) in equation (J3J) since the OZ equations can still be written in the form of equation J5J. The pattern of poles 
is similar to that obtained for hard-sphere mixtures in the Percus-Yevick approximation In figure EH we plot the 
crossover line calculated for size ratio q = 0.5. The packing fractions are given by 77, — pi(Ti, i — s.b, in one dimension; 
Oi is the length of the rod of species i. Crossover occurs from the pole iri to the pole TT2 in a similar fashion to that 
described in figures ^ an d [3 for hard-spheres. Moreover the shape of the crossover line is reminiscent of that in figure 
n Of course the line can extend to r\h — 1 since there is no freezing transition in one-dimension and the maximum 
packing fraction is unity. 

These results demonstrate that crossover occurs in an exactly solvable model. In the vicinity of a crossover point 
we find similar features in the total pair correlation functions as those found in the three dimensional hard-sphere 
mixtures. Recall that for hard-rod mixtures it is possible to calculate wall-fluid density profiles and the solvation 
force exactly (to arbitrary numerical precision) which allows one to ascertain precisely the manifestation of structural 
crossover in these quantities 26] . 



B. Two-Dimensional Mixtures of Hard-Disks 



Theories of liquids are, generally speaking, more difficult to construct and implement in two dimensions than in 
one or three dimensions. This is certainly the case for the fundamental measure DFT that we have used throughout 
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FIG. 9: Logarithm of density profiles of component 6 in a polydisperse hard-sphere mixture, with the particle size distributions 
shown in figure|H| adsorbed at a planar hard- wall. The component b considered here refers to that with a — a™ ax , the diameter 
corresponding to the maximum packing fraction in the distribution of the bigger spheres. For all three curves ^ rjb.i = 0.1. 
In curve (a) the sum of packing fractions of the smaller components ^2 i Vs,i = 0.1 and the wavelength of oscillations is about 
a™ 1 . By contrast, in curve (c), where J^.^s.i = 0.2, the wavelength is about a™ ax /2. In the intermediate case, curve (b), 
with J^ 7 ?*.* = 0.15 we find interference effects. Results should be compared with these in figure[3]for a binary mixture with 
q = 0.5.' 

this paper. Whereas in d = 1 the exact functional for a hard-rod mixture is known and in d = 3 very accurate 

versions of FMT for hard-sphere mixtures [8|, |3, |27J are available, some pre-requisites for FMT-like functionals are not 
fulfilled in d = 2 [2S|. The starting point for FMT is the deconvolution of the Mayer-/ function that describes the 
interaction between two particles, into geometrical weight functions. In even dimensions, however, this deconvolution 
can be achieved only with an infinite number of weight functions. Hence, in practice, the deconvolution of the Mayer-/ 
function in d — 2 can be carried out only approximately (25|. For this reason the FMT results can be less accurate 
than in d = 3, although there have been few careful tests of the theory for mixtures. 

On the other hand, the study of correlation functions in d = 2 is very interesting from an experimental point 
of view. Recent improvements in techniques of video microscopy enable accurate measurement of pair correlation 
functions of colloidal suspensions to be performed in real space pol l2ll |22| . In these experiments, the colloids are 
restricted to move in a plane by the use of light fields, i.e. the system is effectively two dimensional. Experiments 
with one-component colloidal suspensions have demonstrated that the pair correlation function can be measured very 
accurately at small and intermediate separations and, in some cases, even at fairly long range. An extension of such 
experiments to binary colloidal mixtures should be straightforward. 

In order to study the crossover behaviour of correlation functions in d = 2 we performed DFT calculations and 
standard MC simulations in the canonical ensemble for a binary mixture of hard disks with a size ratio of q = 0.5. 
Within DFT we used the test-particle approach, fixing a single big disk at the origin and calculating the density 
profiles pi(r) of both components, i = s,b, in the external field exerted by the fixed disk. From the density profiles, 
we determined the correlation functions hu(r) — Pi{r)/pi — 1. This is the same procedure as in paper I, except now 
we employ Rosenfeld's d = 2 'interpolation' form of the free energy functional |25|. Throughout the calculation and 
simulation we kept the packing fraction of the big disks constant at r\b = pbTrcr^/4 = 0.3. We increased the packing 
fraction of the small disks from a very small value, rj s = 0.075, to a high value, rj s = 0.3. The resulting total correlation 
functions hbb(r) obtained from DFT, (a), and MC simulations, (b) are shown in figure ITT1 The overall agreement 
between these two sets of results is encouraging. We find for the lowest value of r) s that the wavelength of the decay of 
hbb(r) is set by <Jb, the diameter of the big disks, and for the highest value of r] s that the wavelength is set by a s , the 
diameter of the small disks. The DFT results for intermediate values of ij s show evidence of interference and indicate 
clearly that crossover occurs near rj s = 0.1975. The MC results show complex interference effects at the intermediate 
values of r] s . Although we did not determine the precise location of the crossover line in d = 2, as we did in d = 1 and 
d = 3, our results demonstrate that structural crossover occurs in d = 2. 

The fact that crossover is reflected in the oscillatory behaviour for intermediate separation and can be seen in MC 
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FIG. 10: Crossover line for a one-dimensional binary hard-rod mixture with size ratio q = 0.5. On each side of the line 
the asymptotic behaviour is dominated by different leading-order poles 7ri and 7T2. When 7Ti dominates the wavelength of 
oscillations in hij(r) is ~ at whereas when 7T2 dominates it is ~ 0.5<T(,. These results were obtained using the exact result for 
5j?(r) to determine the poles. The dashed line corresponds to total packing fraction rj = r]b + T) s = 1. 



data, which contain some statistical noise, suggests that video microscopy in d = 2 should be well suited to study 
various aspects of the decay of correlation functions in colloidal mixtures. 

VI. DISCUSSION 

We have investigated the nature of structural crossover in various properties of binary hard-sphere mixtures. 
Whereas in paper I we focused only on the decay of pair correlation functions, hij(r), in homogeneous (bulk) mix- 
tures, here we consider manifestations of crossover for the one-body density profiles, Pi{z), of the mixture adsorbed at 
planar walls for the solvation force, f s (L), and for the depletion potential, W(r), of the same mixture. In all of these 
quantities we find a clear signature of structural crossover, i.e. the wavelength of the oscillatory decay changes from 
being approximately the diameter of the big species to approximately that of the smaller species at some crossover 
point in the (^b^s) phase diagram. We argue that the precise crossover point must be the same as that obtained from 
analysing the leading order poles of hij(k), as described in I. Moreover, we show that leading order asymptotics, as 
described by a single pole approximation such as equations J7J , (|13fl and (|16fl provides a remarkably good description 
of oscillatory structure at intermediate range as well as at longest range — provided one is not too close to the 
crossover point, where a two-pole approximation such as © is required. This conclusion, which is equivalent to that 
ascertained in I for hij(r), is important for practical purposes. Experiments or computer simulations are necessarily 
restricted to small or intermediate length scales — one cannot probe the ultimate decay of structural properties. For 
example, if one were seeking evidence for structural crossover in measurements of the solvation force for an asymmetric 
binary mixture confined between the two mica plates of the surface force apparatus j2^| one would have to probe 
plate separations on the scale of a few particle diameters — as in figure Similar remarks pertain to experimental 
measurements |29| of the depletion potential between two very big colloidal particles immersed in a binary mixture 
of big and small (hard-sphere) colloids. Crossover is already apparent at intermediate length scales — see figure |SJ 

It is important to enquire whether the results we obtain here for the idealised system of hard spheres are relevant 
to real fluids, which could be atomic mixtures or mixtures of two types of colloidal particles. This issue was addressed 
in I for bulk mixtures and the remarks are equally pertinent here. Assuming that dispersion forces (power-law decay) 
play a minor role, which could be engineered by suitable refractive index matching in colloidal fluids, structural 
crossover should occur for mixtures where the effective pair potentials are softer than hard-sphe res, provided the 
effective diameters are sufficiently different. Indeed the first observation of such crossover was found |30j in the binary 
Gaussian core mixture where the pair potentials are very soft. Of course, the presence of attractive interactions might 
force the regime where crossover effects can be observed to larger distances. This is the case for the binary mixtures 
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FIG. 11: Logarithm of the total correlation function hbb(r) for a two-dimensional binary mixture of hard disks, with size ratio 
q — 0.5, obtained from DFT (a) and MC simulations (b). rjt = p^a^/A is fixed at 0.3 in all cases. For the packing fraction 
of the small disks rj s = 0.075, the oscillations have a wavelength of about at whereas for r/ s = 0.30 the wavelength is about 
o" s = Ob/2. Interference effects are visible for intermediate values, especially for r/ s = 0.1975 in the DFT results (a), providing 
clear indication of structural crossover. The results are shifted vertically for clarity of display. 

adsorbed at planar walls that exhibit attractive wall-fluid potentials — see figure 

Our results for a simple model of polydisperse hard spheres adsorbed at a planar hard wall warrant further comment. 
We argue, in section llVl that for a size distribution with two broad, weakly overlapping maxima, structural crossover 
remains visible in the density profile of the species b corresponding to the maximum packing fraction of the bigger 
spheres — see figure [3] At first sight this might appear surprising. Consider first a distribution of big hard-sphere 
diameters with a single maximum at a™ ax . Provided the number of diameters is finite, corresponding to a finite 
number of species is, the OZ equations can still be written in the form of equation @ but now generalised to v 
components. Nevertheless, hij(k), with i, j = 1 . . . v, will have poles determined by the zeros of the appropriate v 
component common denominator D(k). There will be a LOP and the quantities ao and a\ will be uniquely determined. 
For a narrow distribution these quantities will be close to those for a one-component fluid with diameter a™ ax at 
the same (total) packing fraction. Similarly for a fluid with a distribution of v small hard-sphere diameters, with a 
narrow peak at <r™ ax , we expect the LOP to be close to that for the one-component fluid with diameter <j™ ax . When 
we consider the polydisperse 2v component mixture, corresponding to the bimodal situation of figure equation (0) 
remains valid and we should expect to find that the LOPs, determined by the zeros of the 1v denominator, are not 
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very different from those in the binary mixture which corresponds to the size ratio a™ ax /<j™ ax . Although v might be 
large (we consider v = 21 in section HV|) . provided it is finite and the ratio a™ ax ja™ ax remains sufficiently small, we 
do expect to find crossover. Of course, as described in section llVl the crossover need not occur at the same point in 
the phase diagram as in the binary mixture. 

We have not attempted to incorporate genuine polydispersity, i.e. a continuous distribution of diameters correspond- 
ing to an infinite number of species, into our theoretical treatment. We leave this for future study, noting that detailed 
DFT and Monte Carlo studies |^, have investigated the density profiles of polydisperse hard-spheres at a planar 
hard-wall. However, these studies employ a parent bulk distribution (top-hat or Schulz) that has a single maximum 
and are not concerned with the intermediate range oscillatory behaviour of the profiles. In experimental samples, 
typical degrees of size-polydispersity are around 5%, which is less than the degree of polydispersity considered here. 
Therefore we are confident that the structural crossover can be observed in experiments with colloidal mixtures. 

Finally we return to our results for one-dimensional hard-rods (subsection IV All . It is gratifying that structural 
crossover occurs in an exact statistical mechanical treatment and that the overall shape of the crossover line follows 
that calculated for the binary mixture of hard spheres. This attests to the generic nature of the crossover mechanism. 
In both one and three dimensional fluids with q = 0.5 the line begins at the origin of the (Tjb,T} s ) plane. (Note that our 
results in figure 5 of I did not extend to very low packing fractions.) For a very low concentration of small particles 
crossover occurs rapidly on adding a small amount of big particles. The form of the crossover line in the limit 775 — > 1 
for the hard-rod mixture in figure ITUI is somewhat more curious. If the crossover line exists in this limit it must, of 
course, terminate at the point r]b = 1, Tj 8 = since the total packing fraction cannot exceed unity. Nevertheless, it is 
striking that for a very large concentration of big rods, say 77^ = 0.97, adding a very small amount of small particles 
so that the total packing fraction is about 0.98, already leads to crossover, i.e. the wavelength of the longest range of 
oscillations in the correlation functions is then determined by the length of the small rods which are present in only 
very small concentrations. We should note, however, that as r\h — > 1 several poles begin to contribute to hij(r) at 
long -range; the imaginary parts of the higher-order poles become comparable with ao, the imaginary part of the LOP 

M _ 

Michael Fisher and Ben Widom 1] pioneered investigations of how the character of the asymptotic decay of the 
pair correlation function should change as one moves from one region of the phase diagram to another. Their study 
is relevant to fluids in which there are both repulsive and attractive contributions to the interparticle potential and 
the existence of the FW line is intimately linked to the presence of liquid-gas, or in the case of mixtures, fluid-fluid 
phase separation. The structural crossover we describe here and in I arises in fluid mixtures where the interparticle 
potentials are purely repulsive; what is required is that the two species are sufficiently different in size. As mentioned 
in the introduction, FW based their study on exactly solvable one-dimensional models. For reasons that will be 
abundantly clear to the reader, we are grateful that they did not consider mixtures! 

R.E. thanks S. Dietrich and his colleagues for their hospitality during his visits to the MPI in Stuttgart. 
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